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Abstract: We present a fully analytic calculation of the Higgs boson transverse mo- 
mentum and rapidity distributions, for nonzero Higgs at next-to-leading order in the 
infinite-top-mass approximation. We separate the cross section into a part that contains 
the dominant soft, virtual, collinear, and small-p^-enhanced contributions, and the remain- 
der, which is organized by the contributions due to different parton helicities. We use this 
cross section to investigate analytically the small-p_L limit and compare with the expecta- 
tion from the resummation of large logarithms of the type hymn /pl- We also compute 
numerically the cross section at moderate px. where a fixed-order calculation is reliable. We 
find a K-ia,ciox that varies from ~ 1.6 — 1.8, and a reduction in the scale dependence, as 
compared to leading order. Our analysis suggests that the contribution of current parton 
distributions to the total uncertainty on this cross section at the LHC is probably less than 
that due to uncalculated higher orders. 



Keywords: 



QCD, Higgs Physics, Hadronic Colliders. 



Contents 



|l]. Introduction ^ 

|2|. Higgs p± Spectrum: General Formulae and Leading Order expressions ^ 

H. 0{as) Corrections g 

|3.1| Virtual Contributions ^ 

3.2 Altarelli-Parisi Contributions P 

3.3 Real Contributions 
|3.4| Total Corrections ^ 

^ Numerical Results and Analysis 12 

H]. Small-px limit [l4 

|6|. Conclusions 17 

|A| . Nonsingular Real Contributions 18 

[All 99 ^ H + X terms |l8 



A]| gq^ H + X terms || 



A.3| QiQi ^ i7 + X terms 

A.4 Qj^j ^ i7 + X terms ^ 
ATS] gi^i ^ + X terms |§ 

g. Phase space for integrals over momentum fractions 

|Q. Details of the small-p^ limit 27 



1. Introduction 

Our current understanding of particle physics depends crucially on the breaking of the 
electroweak symmetry to give masses to the and Z bosons, as well as to all of the 
matter particles. Yet, despite decades of extremely precise and successful predictions, 
the exact mechanism by which this symmetry breaking occurs is not known. The simplest 
model involves a single weak SU{2) scalar doublet which has a nonzero vacuum expectation 
value. After rewriting in terms of the physical states, this "standard model" leaves behind 
a single neutral scalar, the Higgs boson, as its signature Q. Furthermore, models with 
extended Higgs sectors, including the minimal supersymmetric standard model (MSSM), 
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often have a relatively light scalar with properties similar to the standard model Higgs 
boson. Therefore, the Higgs boson is at the top of any list of new particles to be found. 

The direct search for the standard model Higgs boson in the e~^e~ — > HZ channel at 
LEP2 has put a lower bound on its mass of 114.1 GeV [|[. Moreover, there are hints of a 
signal in the data just above this bound. Meanwhile, precision electroweak measurements 
give an upper limit of mn ^ 196 — 230 GeV at the 95% confidence level Therefore, 
if a standard model Higgs boson exists, its allowed mass range is not large. Run II of the 
Tevatron can exclude a standard model Higgs boson over much of this mass range, up to 
about 180 GeV, assuming 15 fb~^ per experiment. However, a definitive 5a discovery is 
difficult to obtain at this luminosity for a mass much beyond the LEP2 limit. The Large 
Hadron Collider (LHC) at CERN will be needed to certify a Higgs discovery, and to pin 
down its mass and couplings. 

The dominant production mechanism for the Higgs boson at the LHC is the gluon-gluon 
fusion process. This process occurs at leading order (LO), O(a^), through a heavy quark 
loop. As is typical in Quantum Chromodynamic (QCD) processes initiated by gluons, the 
radiative corrections are quite large. The next-to-leading order (NLO) corrections have 
been computed, including the full top-mass dependence |^, and indeed, the ET-factor is 
~ 1.4 — 2.2, depending on the Higgs mass and the scale choice. The computation with full 
mj-dependence requires the calculation of two-loop diagrams and is quite complex. Luckily, 
it simplifies greatly in the limit of large top quark mass {nit — > oo). In this limit, one can 
integrate out the top quark loop, leaving behind an effective gauge-invariant Hgg vertex Q. 
The Higgs boson production cross section has been calculated at NLO in this limit in 
Refs. [^. It gives an excellent approximation to the full mt-dependent NLO cross section 
for niH ^ 2mt. Furthermore, the JC-factor calculated in the effective theory gives a good 
approximation to the full mj-dependent if -factor, even for larger Higgs masses. Attempts 
have been made to estimate the NNLO corrections using soft-gluon approximations |Q and 
resummations Q. Recently the full NNLO cross section in the large-mj limit has been 
computed |l^. Although the K factor is larger still at NNLO, the increase is not as severe 
as the NLO enhancement, and the perturbation series seems to be well-behaved. 

In addition to increasing the total cross section, the QCD radiation can have a large 
effect on the kinematic distributions of the Higgs boson. Most notably, the transverse 
momentum, p±, of the Higgs boson is exactly zero at LO, but is typically nonzero at 
higher orders, due to additional radiated partons. In fact, this additional QCD radiation 
has led some to consider searching for the Higgs boson in association with a tagged jet at the 
LHC [ pT|] . Regardless of whether additional jets are tagged in the Higgs events, it is useful 
to understand the transverse momentum distributions of both the Higgs signal and the 
background. The transverse momentum spectrum of the Higgs boson has been calculated 
at LO*, both with the full top quark mass dependence and in the large- mt limit |12]. It 
was seen that the large-mj limit is a good approximation to this distribution if niH ^ 2mj 
and p± < mt- Recently, it was shown [^3| that these conditions are also sufficient to use 
the large-mj limit for Higgs +2 jet production; i.e, the transverse momenta of the Higgs 

*Note that the LO contribution to the p± spectrum at non-zero p± is actually down by Us compared to 
the LO cross section; it contributes to the cross section at NLO. 
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boson and the jets must be less than ~ mt, but other invariants such as the total partonic 
center-of-mass energy can still be large. 

In this work we calculate the NLO corrections in the large-mj limit to the Higgs boson 
p± and rapidity spectrum, da/dp'^ /dyn, at the LHC. Initial results of our calculation, which 
included the purely gluonic contributions, were reported in Ref. |14|. The NLO corrections 
in this limit have been calculated previously in Ref. [^], using a Monte Carlo integration 
to do the phase space integrations, after cancellation of the infrared singularities. More 
recently, a full analytic calculation of the p± spectrum was presented in Ref. ||l^. Our 
calculation is similar to that in |jl^ in that it is completely analytic. However, by treating 
the contributions of different helicity configurations separately, we are able to report the 
full formulae for the differential cross section in this paper, in a relatively tractable form. 
Another difference is that we have used the +function technique to deal with soft singular- 
ities, whereas an artificial parameter was introduced in |16| to separate the soft and hard 
radiative contributions. Our technique allows more of the universal structure of the NLO 
corrections to be apparent. Finally, since much of the total Higgs cross section occurs at 
not-too-large p±, we investigate the low p± limit of our result. 

At small Higgs transverse momenta the perturbation series for the p± spectrum be- 
comes unstable, containing terms like (a^/p^)a^ ln'"(m|^/p^), with the leading logarithm 
occurring for m = 2n — l. This logarithmic series has been resummed, using the techniques 
of Collins, Soper, and Sterman at various levels of approximation The NLO 

differential cross section at small p± contains the fixed order terms in the logarithmic ex- 
pansion corresponding to n = 2, m = 3, 2, 1, 0. In particular, the so-called B^"^^ coefficient 
in the logarithmic series occurs in our calculation at n = 2, m = 0. This coefficient has 
recently been derived |2^ using the universality of the real emission contributions, com- 
bined with knowledge of the virtual correction amplitudes in the soft and collinear limits. 
Using our analytic expressions for the NLO p± spectrum, we have verified this coefficient 
by direct calculation. The analytic comparison of our cross section in the small-pj_ limit 
with that expected from the resummation formulae is a very stringent check on our results. 

The remainder of this paper is organized as follows. In section ^ we set up the cal- 
culation by defining some general formulae and giving the Born level expressions for the 
differential p± and rapidity spectrum. In section |3| we obtain the 0(a<j) corrections to 
this by combining the virtual one-loop (in the effective large-mj theory) amplitudes with 
the tree-level real radiative corrections. Although both of these contributions are infrared 
divergent, these divergences cancel after they are added together, using MS parton density 
functions, defined at NLO. In the main text we give formulae for the largest contribution 
to the distribution, which contains all terms having singular behavior as one of the real 
QCD partons becomes soft or collinear, as well as most of the terms that are leading at 
small-pj_ (all but the m = terms from the previous paragraph). For lack of a better 
word, we label these contributions the "singular" contributions. The remaining "nonsin- 
gular" contributions are given in appendix In section Q we give some numerical results 
and analysis obtained from our calculation, showing some representative distributions at 
the LHC. We also comment on the numerical comparison of our results to previous calcu- 
lations. In section |5| we consider the small-p^ limit of our result, and compute directly the 



-3- 



B^'^^ coefficient. Finally, in section ^ we give our conclusions. 



2. Higgs p± Spectrum: General Formulae and Leading Order expressions 

In the large-mj limit the top quark can be removed from the full theory, leaving a residual 
Higgs-gluon coupling term in the lagrangian of the effective theory: 



1 -— 1 + -^A 

Svr u V 47r 



(2.1) 



The finite 0{as) correction to this effective Hgg operator, which is necessary to the order 
we are working, has been calculated to be 



57V, - 3Cp 



11 



(2.2) 



where iVc = 3 and Cp = {N^ - l)/(2iVc) = 4/3. The ©(a^) correction to the effective 
operator has also been calculated ||l^, ^. 

We consider the inclusive scattering of two hadrons hi and /12 into a Higgs boson, 
/11/12 H + X. The differential cross section for this process, via the gluon-fusion produc- 
tion mechanism in perturbative QCD, for a Higgs boson of transverse momentum p± and 
rapidity yn, can be written 



da 



dp]_dyH 



dXa dXb fi/h^ (Xa, HF)fj/h2 i^b, /"f) 



da. 



dpidyn 



(2.3) 



where i and j label the massless partons which scatter from the hadrons hi and respec- 
tively (i = g, qf,qf, with nj flavors of light quarks) in the partonic subprocess ij H + X. 
The parton densities fi/h{x,iiF) are defined in the MS factorization scheme at scale fip. 
The partonic subprocess cross section can be expanded as a power series in the strong 
coupling, asinR), as 



ddij 
dp^dyu 



<7o 



as(/i 



«V(i) 



27r 



2tt 



R) 



G. 



(2) 



+ 



(2.4) 



with 



64 V Sttv J 



(2.5) 



where the Higgs vacuum expectation value v is related to the Fermi constant, Gp = 
1.16639 X 10~^ GeV~^, by v~'^ = \/2Gf- The quantities G^j ^ are functions of the renor- 
malization scale fiR, the factorization scale ^p, the Higgs mass mn, and the partonic 
Mandelstam invariants, defined by 



S = {Pa+Pb? = {PH + Qf 
i= {Pa-Q? = [Ph-PH? 
U = (Pb- Qf = [Pa - Puf 



SXuXb 

m\ — \Gxbm±e^'' 



(2.6) 
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where pa = XaPi and pi, = XbP2 are the initial-state parton momenta, Q is the momentum 
of the final-state partons which balance the Higgs boson, and = m'jj + p^- These 
invariants satisfy 

g2 = s + i + u-rriH . (2.7) 

At leading order the contributions of the different subprocesses (at nonzero p_l) are 
given by 



with 

9gg = N, 



mf^ + s'^ + + 
uts 



9a, = Cf ( ) (2.9) 



Qqq — 2Cp 



and Qqg obtained from ggg via i ^ u. Note that the quark and antiquark in gqq must be of 
the same flavor at leading order. 

3. 0{as) Corrections 

The 0{as) corrections come from three sources: the virtual corrections (V) to Higgs-plus- 
one-parton production arising from the interference between the Born and the one-loop 
amplitudes; the real corrections (R) from Higgs-plus-two-parton production at the Born- 
level; and the Altarelli-Parisi corrections (AP) arising from the definition of the MS parton 
densities at NLO. 

Although the total 0{as) correction to the inclusive Higgs production cross section 
is finite and well-defined in four dimensions, the individual V, R, and AP contributions 
are not and must be calculated using some regularization procedure. As is standard in 
perturbative QCD we use conventional dimensional regularization (CDR) with dimension 
d = 4 — 2e, so that the collinear and soft singularities appear as e-poles in the separate 
contributions. The total 0{as) correction is 

Gg) = lim(Gf )(e) + Glf)(e) + G;f^)(e)) . (3.1) 

3.1 Virtual Contributions 

The virtual one- loop expressions in the effective large-mt theory were calculated in Ref. . 
Although they were calculated for helicity amplitudes in the t'Hooft-Veltman (HV) and 
the four-dimensional-helicity schemes [El], they can be utilized in CDR by using the results 



of Catani, Seymour and Trocsanyi [24|. In that paper it was shown that the singularity 
structure at NLO in the HV scheme and the CDR scheme are identical, when properly 
defined, because the number of internal gluon helicities are taken to be the same; i.e., in 
d = 4 — 2e dimensions. This implies that we can reconstruct the CDR virtual contributions 



- 5 - 



by first constructing the virtual corrections to the cross section in the HV scheme, using 



the amplitudes from Ref. |21] with 6r = 1, and then replacing the 4-dimensional Born- 
level expressions gij ( |2.9[ ), which multiply all virtual singularities in the HV scheme, by the 
following d-dimensional generalizations: 



1 - e)(m|^ + + r + - Aemjjiits 



uts 



1 - e 



(3.2) 



gqq{e) =a,2Cp ' 



where 



2 \ e 



r(i -e)\ pI J ■ 



(3.3) 



This accounts for the correct number of helicities of the external particles in the CDR 
scheme, as well as the d-dimensional final-state phase space. 

3.2 Altarelli-Parisi Contributions 

The Altarelli-Parisi subtraction terms arise from the renormalization of the parton density 
functions at NLO. In the MS scheme they can be obtained from 



er(l^- e) (^^) ^^Pk^{z)9kJ{zXa,x^,■,e)6(i+{Q^ -i)z'^ 

+Pkj{z)gik{xa,zxh]e)5{u + {Q^ -u)z^ > , (3.4) 



where we have made explicit the functional dependence on Xa and Xh of the gij{e), given 
in (|3.2|). The splitting functions are 





= 




= Cf 


Pggi.^) 


= Cf 


Pggi^) 


4[ 



l + z^ + jl-z) 
(1 - z)+z 



l + (l-2)2n 



+ (3oS{l - z) 



(3.5) 



[z' + (1 - zf 



where (5q = (11 Ac — 2n/)/6 and is the number of light quarks. The -l-functions are 
defined by 



C dzF{z){g{z))^ = C dz[F{z)Q{z - x) - F{l)]g{z) 

Jx Jo 



(3.6) 
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where in eqs. (|3.5| ) we have g{z) = (l — z) ^. The integrals over z in ( ^.41 ) can be performed 
using the 5- functions, giving 



where 



Pki{za)9kj,a{^a; e) + Pkj{zb)gik,b{zb; e) 



-t 



(3.7) 



—u 

= W^u ' ^^-^^ 

and we have defined gij^a[za-,e) = gij{zaXa, X),; e) and gij^b{zb]() = 9ij{xa, ZbXt; e). Note 
that gij^a{za;e) can be obtained easily from the formulae of (^) by the replacements 
s — > ZaS, t ^ t, u ^ ZaSp\/t. The expression gij^bizb]^) can be obtained analogously. For 
later use, we define gij^aiza) and gij^bi^b) as the e — > limit of these. 

3.3 Real Contributions 

The amplitudes for the relevant Higgs-plus-two-parton-production processes at the Born 



level in the large- mt limit were calculated in Refs. |25, |2^, using the helicity amplitude 



method^. Since we are using the CDR scheme at NLO, we must supplement the resulting 
4-dimensional squared-matrix elements, \M{d = 4)p (averaged over initial polarizations 
and colors and summed over final ones), with the correct e-dependence to obtain the 
squared-matrix elements in d = 4 — 2e dimensions, \M{d = 4 — 2e)p; i.e., we must add the 
quantity |A4((i = 4 — 2e)p — |7W((i = 4)p. In practice this was done by investigating the 
collinear regions of the phase-space, where correction terms, finite as e — > 0, were obtained 



in accordance with Ref. |24|. Most of these correction terms combine so that they may 
be identified with 0{e) parts of the splitting functions, or universal terms multiplied by 
5{Q'^). The single remaining form of correction term, which is not so easily identified, turns 
out to be crucial for obtaining the correct small-p^ limit. Thus, we are confident that our 
approach gives the correct CDR cross section at NLO. 

The real contributions from initial-states i and j can be written 




(^)"^->-^fl^}'<-) 

where 6 and (j) are the polar angles of the recoiling partons in their center-of-mass frame, 
and in addition to averages and sums over colors and polarizations, we also sum over 

t There are several typographical errors in the formulae for the squared matrix elements in the appendix 
of Ref. j26|, which we list here. The occurrences of in the denominator of eq. (A12) should be replaced 
by 524. The expressions for nis and 1123 in eq. (A16) should be replaced by nis = — S'24ni2(l ^ 4) and 
?i23 = ?^l3(3 ^ 4). The interference term that occurs in the scattering of identical quark pairs has the 
wrong sign on the right hand side of eq. (A20). 
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possible final-state flavors. By performing many partial fractionings on the squared- matrix 
elements, we can reduce the integrals to two types [27]: 



(m,n) 
kl 



n 



1 

2^ JO 

(m,n) 1 



akl 



de (sine)i-2^(sin<^)-2^ 



(3.10) 
(3.11) 



where Sak = (Pa —PkY ^'^^ ^abk = ipa + Pb — Pk)'^ with a, b labeling the initial-state parton 
momenta and k, I labeling the final-state parton momenta. A general solution to integrals 
of type ( |3.10| ) for arbitrary e was given in Ref. |28]. Solutions to the integrals of type 



( 3.11 ) were given as a power series in e in Ref. |27]. We have independently checked these 



integrals, and supplemented some of them to retain terms of 0{e) not given in |27|. 

The integral of the real-emission partonic cross section over the momentum fractions 
Xa and Xb in eq. (2^) has infrared singularities as 0. To handle these, we used the 

following identity: 



-l—ne 



ne 



-i/za)' 



(1 - Zah 



ne 



, (3.12) 



and the equivalent equation obtained by a ^ 6 and t ^ u. 

To organize the final expressions, we write the real contributions as a sum of "singular" 
and "nonsingular" terms: 



(2R,ns) 



(3.13) 



where the "nonsingular" terms are finite as e ^ and as 0. We also choose to 

define the separation so that they have no naive singularities as — > as well, prior to 
the integration over Xa and Xb in Eq. ( |2.3D ; this ensures that the "nonsingular" terms are 
not enhanced by any power of ln(m|^/p^) in the small-pj_ limit. (They will have l/pj_ 
singularities after the integration over Xa and Xb , but without logarithmic enhancement.) 
Explicit formulae for the "nonsingular" terms are given in appendix 

3.4 Total Corrections 

After combining the three contributions, the e-poles cancel, and we can safely set the 
dimensional regularization parameter e to zero. We can then write 



G 



(2s) 

where the "singular" terms, G-- , contain the sum of the virtual, the Altarelli-Parisi, and 



(2) 



G 



{2R,ns) 



(3.14) 



(2s) 

the "singular" real contributions. In this section we give formulae for the G • , which 



contain by far the dominant contribution to the cross section. In the small p± limit they 



contain all of the contributions enhanced by (ag/p^)a^ ln'"(m|^/p^), with 



m 



3,2,1. 



The "nonsingular" terms, G- , which we give in appendix contain only subleading 



m 



contributions in this limit 
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In the following equations we use the quantity Q\ = Q"^ +p\ and the definitions 



and 



Pggi^) = (1 - z)Pgg 

Vqq{z) = {l-z)Pqq 



Nr 



(3.15) 



= C^(l-z) 
Clq{z) = CFZ 
C|,(z) = ^(l-z). 

Note that the Cfj are minus the 0{e) parts of the splitting functions. 
For the gluon-gluon singular term we obtain 



99 



d{Q^)}^{A + S + N,U)ggg 

+{Nc - uf)^ [{m%/s) + {m%/i) + (m|,/u) + ml] 



M 



1 



-Pgg{Za) In + Pgg{Za) (^^yr^ 



A*>^a . , , /Inl 



ggg,a{Za) 



+ 



2n/P,,(2;„)lngt + 2n^C^^(^„) 
lnl-z,\ In ' 



9q9,a{Za) 



(1 - ^a) + 



(3.16) 



8^ + |4 + g8 + ^4 + £4) + ^^^^ (^8^ + 5^ + + (^i/^b)^ + (f/^^)^) ' 



silt 



1 ^ /?q ^ {ml^s^ + ZaZb {{u/zj,Y + (f/^a)^) 



-tj \l- Za 



silt 



+ 



+N, 



where 



{i, a) ^ (u, b) I 

(ml, + s^ + Q^ + {ii/z^Y + ii/zg)^) (Q^ + Qi) 
s^Q^Ql 

2m% {{ml - t)4 + {ml - ii)* + + t^) 
silt {ml — ii){ml — t) 



2 \ -t Vis 9 ^ 



1 1 PI 



^ > (3-17) 



(3.18) 
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and 



1,2-"" 

2 -t 3 

— In -At- In — ^ In -At- In „ 



U —t —u 

— In — TT- In 



2 2 



+ In2!!^ + ln2 



m 



H 



4 



+ ln^ 



H 



mjj — u 



+ 2Li2 



f s — mjj 



+ 2Li2 



— t 



+ 2 Lis 



(3.19) 



(3.20) 



Here, Lis is the dilogarithm function. 
The gluon-quark singular term is 



/7(2s) 



= (5(Q2) } (A + iVeFl + CfFs + 1^3) ggq 



+{Nc - Cf)Cf 

1 



-t 



+ P + — ■um|^ 
— n 



o / X 1 M^^a , / / In 1 - 



+ 



+ 



+ 



1 



-U 



-u 



,2 



5'g?,o('^a) 
In 1 — 



9gq,b{Zb) 



In 1 — Zn 



9gg,b{^b) 
In {QlzJi-i))\ 



xN.Cf 



^- Za J + (1 - ^;a)+ J 

■sH-si^ + qH+qH^ 
sut 

ZaZb {-S^ii/Za) " s{i/Zaf " ^'^(ii/^b) - Q^{u/Zbf) 



silt 



Zb 



3^2 /*' + ^' 



(-53(i7zJ - S(i7z„)3 - Q^iu/Zh) - Q2(^i/^,)3) (g2 + g^^ 



2m|^ ((m^^ - t)2 + P 
su{m'\j — u) 



^In^ 

pI Ql 



(3.21) 
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with 



2 -t 2 -u 2 -t 

^ s ^ —i , I , —u 
+ In — 2~ — 2 — T — T 



In — TT- In — ^ 



my my 



+ 2Lj^^Uln^^^ + vr^ 
\my — u) my — u 



Vo = In^ ""^ ' 



+ ln' 



m^jj — t 



2 In At- In ^ 



H 



+ 2Li.(i^) +2Li,( 



mj^ — t 



7 _ 27r2 
2 3 



^3 = /3o(21n^+ln^ 
— u -i 



67,, 10 



(3.22) 



(3.23) 
(3.24) 



For the quark-antiquark terms we must now distinguish between the scattering of 
identical or nonidentical flavors. The quark-antiquark singular term (same flavor) is 

G^l = (5(Q2) J (A + iVe^i + CfW2 + W^) 



+{Nc - Cf) 2Cl 
1 



-t 



P + + — sm\ 



-PM In ^ + P,,{Za) [fr^ ) ^ + CI,{Za) 



+ 



Inl - Za 



-^..(^a)ln^+q,(^a) 



+ 



9gq,a{Za) 

Za \ ( f Inl- Za\ In {Q'iza/i-i)) \ 



ij \ \ 'i-- Za J + (1 - Za)+ J 

P + u2 + (^i/Za^ + {U/Zbf 



X {2Cf - N,) Cl 

1 



-tj \l- Za 



9qq,a{Za) 



+ 



(t, a) ^ (u, 6) 



+2C| 



{s-Q'')^ + {u + i-2Q^y 



1 , Pi 



(3.25) 



with 



Wi = In — ^ hi — = hi — ^ hi 



m^ m\j 



2 2 

my 



1" ^ 1^ ZX 
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+ 2Li2 



s — m 



H 



+ In 



2 rriH 



2 -t 



2 — 



1 1 s 
m — ; + m 

—t —u 



u 



2 In In 



In' 



m 



H 



my — u 



In^ 



m 



H 



mjj — t 



+ 2Li2 



ni 



H 



mjj — u 



+ 2Li2 



m 



H 



mjj — t 



H 



7 + 2tt' 



Finally, the quark-antiquark singular term (different ffavors) is 
.(2s) I / 1 



G. 



-t 



Q 



9gq,a\Za) 



{t, a) {u, b) 



+2Cl 



{s-Q^f+{u + i-2Q'f' 



1 1 pI 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



The singular terms for the quark-quark scattering, both same flavor or different flavor 
quarks, are identical to ( 3.29| ). 



4. Numerical Results and Analysis 

In this section we present some numerical results of our calculation. A convenient choice 
of variables that we use for the integrals over Xa and X}, is given in appendix ^ We choose 
to make all plots for proton-proton collisions at the LHC center-of-mass energy of ^/s = 14 
TeV, and for a Higgs mass of mn = 120 GeV. This value of the Higgs mass is not far above 
the direct search limits found by LEP2. Unless otherwise indicated, we use as a standard 
choice of renormalization and factorization scales, = hf = m±. For parton density 
functions (PDFs) we use the CTEQ5M1 distributions for all NLO plots and the CTEQ5L 
distributions for all LO plots |2^, with the corresponding values of as{mz) = 0.118 and 
as{mz) = 0.127, respectively (with as running at corresponding order). We discuss the 
consequences of this choice and the error due to PDFs in general below. 

Before presenting our results we must first understand the range of Higgs p± for which 
they are applicable. In Fig. |l| we plot a comparison of the p± spectrum for yn = computed 
at LO, with the exact rrit dependence (dashes), and in the mt = oo effective theory (solid). 
For p± < 200 GeV the effective theory underestimates the cross section by about 5 percent, 
which, as we shall see, is much less than the NLO corrections. Above 200 GeV, the 
exact calculation and the effective theory calculation begin to diverge, with the effective 
theory calculation overestimating the cross section. This behavior is similar for other Higgs 
boson masses as long as tuh ^ 250 GeV, where the effective theory underestimates the 
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cross section by about 25 percent in the low p± region. Thus, we assume that our NLO 
calculation will be good for p± < 200 GeV, and is applicable for the Higgs masses that are 
favored by the electroweak data. Prom here on, all of our plots will be calculated in the 
rrit = oo effective theory. 

The effect of NLO corrections on the Higgs p± and rapidity distributions have been 
studied numerically before in Refs. and [^]. Since these two calculations appear to 
be in agreement^, we have only checked our calculation directly against that of Ref. |16|. 
Using MRST99 PDFs, and rescaling the effective Hgg coupling to be consistent with their 
definition, we have been able to reproduce Figs. 3 and 4 from with excellent agreement. 

In Fig. ^ we plot the Higgs p± spectrum for three different values of the Higgs rapidity 
yn- These curves show the standard features of this distribution: the fall off with higher 
p± and larger yn , and the steep plummet of the curve at very small pj_ <^ rriH ■ This last 
feature is due to the large logarithms \D.{mH /p±)^ which give a large negative contribution 
at NLO at small p^. Although the fixed order perturbation theory is unstable in this 
region, below about p_\_ = 30 GeV, we continue the plot here because we are interested 
in the connection between the fixed-order cross section and the resummation of the large 
logarithms. This will be discussed in more detail in section ^. 

In Fig. m we plot the separate contribution of the different initial-state partons to 
the p^ spectrum at y// = 0. At small p_\_ the cross section is strongly dominated by the 
gluon-gluon component, whereas at large ~ 200 GeV the (anti) quark- gluon contribution 
makes up about 37% of the total cross section. The (anti)quark-(anti)quark contribution, 
summed over all fiavors, is very small and negative over most of the range of this plot. Its 
magnitude is less than one percent of the total for px_ > 30 GeV, where the fixed order 
perturbation theory is reliable. 

In Fig. 1^ we plot a comparison of the LO distribution with the NLO distribution for 
yu = 0, while varying the renormalization and fragmentation scales together (/i/j = fip) 
between 0.5m± and 2m±. This exhibits the two main features of the NLO cross section, 
relative to LO: it is substantially larger, and the uncertainty due to scale variation is 
smaller. We illustrate these two features in more detail in Figs. ^ and ^. In Fig. |5| we plot 
the i^T-factor, defined as the ratio of the differential cross section calculated at NLO and 
LO; that is, 

^ ^ {d^/dpl/dyH)^^o _ ^^^^^ 
{da/dp\/dyH)^Q 

In this plot we have set the scales to /x/j = = rn_\_. We also reiterate that the LO 
cross section in this ratio has been calculated using the LO PDFs CTEQ5L. For > 50 
GeV, the ii'-factor is relatively constant, rising slowly from 1.7 to more than 1.8. This is 
comparable to the NLO ii'-factor for the total cross section, for this value of the Higgs 
mass, mH = 120 GeV. For p± < 50 GeV the K-factor begins to drop more dramatically, 
less than 1.6 at = 30 GeV, where the large logarithms start to become important. 



'''Although the version of Ref. jl6j pubhshed in Nucl. Phys. B 634, 247 (2002) claims some disagreement 
with Ref. the most recently revised version of their preprint, hep-ph/020114v4, which appeared after 
the published version, now agrees with The change was due to a wrong implementation of the MRST 
PDFs in the earlier version of [n6f. 



- 13 - 



Below this, the A'- factor plummets, indicating the inapplicability of this fixed order cross 
section below p± 30 GeV. 

In Fig. 1^ we plot the scale variation of the differential cross section at the particular 
value of transverse momentum pj_ = 50 GeV and rapidity yn = 0. We plot both the 
diagonal variation /i^ = fip = the anti-diagonal variation /i/j = X'^i? = 

m^/x, with the cross section normalized to the value at ^ur = fip = rn±, (x = 1), and 
0.1 < X < 10- For both the diagonal and the anti-diagonal variation, the NLO curves are 
much less sensitive to the scale choice than the LO curves. For instance, for a typical range 
0.5 < X = Ai/j/^T-x = fJ'F/^± < 2, the LO ratio varies from 0.74 to 1.38, whereas the NLO 
ratio only varies from 0.83 to 1.22. In Fig. |6| we also see that the scale variation is least 
when the renormalization and factorization scales are varied together with ^ur = fip. 

In addition to the p± spectrum at fixed yn, we can treat the differential cross section 
as a rapidity distribution for fixed p±. This is plotted in Fig. ^ for p± = 50 and 100 GeV, at 
both LO and NLO. Again we see the steep fall off for large rapidities, due to the restriction 
of the available phase space, and we see the increase of the NLO cross section over LO. 
In Fig. 1^ we plot the ET- factor for p± = 50 as a function of the Higgs rapidity yn- It is 
very flat as a function of yn, except at large rapidities, where it drops slightly, from 1.7 to 
about 1.6 at I y// 1 = 4. 

As we have seen, the NLO corrections give a if-factor of about 1.6 to 1.8 in the range 
of applicability of this fixed order calculation, and this K factor is fairly independent of 
the exact values of pj_ and yn- In addition, the standard variation of scales between 
0.5m_L < fipt = fip < 2m± suggests an uncertainty in this NLO calculation on the order of 
about 20 percent. One may also wonder what the uncertainty on this cross section is due to 
the parton density functions. We have calculated the NLO p± and yH distributions using 
the MRST99 set 1 distribution |^] and found that the difference is less than 4 percent 
over the full range of 30 GeV< p± < 150 GeV and {ynl < 2.5. If we restrict to yH = 0, the 
deviation is never more than 2 percent. This is in contradiction to the results of Ref. |16|, 



where they found a much larger variation, albeit using the older PDFs of MRST98 |]31| , 
CTEQ4 1 32], and GRV98 [33|. Presumably, the more recent global analyses produced less 
variability in their standard PDFs. We expect that the differences resulting from the use of 
the more recent CTEQ6M [34| and MRST2001 [35| will also be small, since neither update 
has a large change in the region of x and gluon distributions relevant to our calculation. 
However, this does not necessarily mean that the uncertainty due to the PDFs is so small. 
The recent analysis on PDF uncertainties in |34| suggests that the expected error on the 
gluon luminosity function in the range of relevant s at the LHC is on the order of 5 percent. 



5. Small-px limit 

At small p± <^ TnH, the p± spectrum becomes unstable at any fixed order, due to large 
logarithms. At LO it diverges to positive infinity as p± — > 0, as seen in Fig. |l], while at 
NLO it diverges to negative infinity, as seen in Fig. ^. These logarithmic corrections can be 
resummed, however, resulting in a well-behaved physical p± spectrum, even for p± <C rrin- 
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In this section we compare our NLO results with what is expected from the resummed 
cross section. 



Following Collins, Soper, and Sterman |17], we can write the resummed Higgs pj_ 
spectrum as an integral over an impact parameter: 



da , mjj ao{asimH)) 



dp\ dyn 



(resummed) 



2s 



bdbJo{bpi)W{h) 



(5.1) 



with 



W{b) = (Cg,{as{bo/b)) o /,) (xO; bo/b) (C,,- (a, (60/6)) o /, ) (x"; bo/b) 



X exp 



H dq^ 



A{as{q))\n^ + B{as{q)) 



(5.2) 



where i and j are implicitly summed over all massless partons {g,qf,qf), and 60 = 2e~'''^, 
with Euler's constant ^e- In eq. (|5.2| ) we have introduced the symbol o to denote convolu- 
tion, defined by 

'1 dz 



In eq. (|5.2| ) the convolutions are evaluated at 



—f{z)g{x/z) . 



(5.3) 



-0 _ mne 



Vh 



^0 



(5.4) 



and the factorization scale in the PDFs is set to bo/b. The parameters A(as), B{as), and 
the convolution functions Cij{as, z) can be expanded as a power series in ag'. 
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The known coefficients and functions in the power series (relevant to Higgs production) 

are [I 



19, 



^(2) = 2Nr 



18 9 ^ ' 



-2/?o 



-25P(,2)+^„(^A + ^Ar,) 



(5.6) 



- 15 - 



A vr^ 
- + -AT, ) 5(1 - z) , 



CS!iz) = C;,{z) = Cfz , 
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where 



^ Nl[- + 3C(3) 



99 



1 2 



(5.7) 



is the coefficient of the (5(1 — z) term in Pgg (z), the NLO splitting function, and C{^) is the 
Riemann zeta function with ^(3) = 1.202057. . .. The comparable coefficients in Drell-Yan 
production have been given in Ref. [37|. All of the coefficients and functions in eq. (p.(]|), 



except 5(2) and C^'g'iz), are universal, in the sense that they only depend on the color 
charges of the initial-state partons involved in the scattering (gg for the present process), 
and not on the final-state that is produced. On the other hand, the expressions for Cg]^ (z) 
and 5(2) given above are specific to Higgs production. Of these coefficients, i?(2) was 
obtained most recently in Ref. |2^ , using the universality of the real emission contributions 
at small p±, combined with knowledge of the virtual correction amplitudes in the soft and 
collinear limits. The universality structure of all of these coefficients was revealed in more 
detail in Ref. where it was shown that the general resummation formula, eq. (^.Ij), 
can be rewritten in such a way to make all of A{as), B{as), and Cij{as,z) universal, by 
extracting the "non-universal parts" into a single new process-dependent parameter H{as)- 
The resummed formula, eq. (^.ID, can be expanded as a power series in ag for < 



p± <^ ran, as in Ref. |39]. At NLO this yields 



dcT 



dp^, dyn 



s pI_ 
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.m=l n=0 
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In^ 

PI 



(5.8) 



where the coefficients, mCn, depend on the coefficients and functions given in eq. (^.6|). 
The formulae for the mCn are given explicitly in appendix^ in eq. (C.2). In particular the 
coefficient depends on B^'^\ We have checked analytically that the small-p_L limit of 
our calculation agrees exactly with eq. (^.^). Thus, we have verified the formula for i?(2) 
for gg ^ H + X given in eq. ( |5.6| ). Some of the details involved in taking the small-p^ 
limit of our cross section are given in appendix 

In Fig. ^ we plot a comparison of our calculation at both LO and NLO versus the 
corresponding small-p_|_ limit formulae, eq. ( f).8\ ). As can be seen from the figure, the small- 
p± limit works very well for p± < 10 GeV, and it is indistinguishable on the plot from the 
exact fixed order calculation for p_l < 5 GeV. 

In Ref. [16 1 the soft-plus- virtual (S+V) approximation, which retains all of the singular 
terms of the cross section as — > 1, was considered as an estimate of the full NLO cross 
section. This approximation had previously been used to estimate the NNLO corrections 
to the total Higgs production cross section In |1€] it was seen that this approximation 
is reasonable for larger p±, where the phase space for non-soft emissions is constricted. 
In this vein we may consider our "singular" contributions from eqs. ( ^.17 ), ( 3.21 ), ( 3.25| ), 
and ( 3.291 ) to be a combined S-l-V and small-p^ approximation, as it contains all of the 
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dominant soft-plus- virtual terms, and it contains all of the leading contributions at small 
p± except the 2C0 terms. In addition, we note that it also contains the contribution of all 
collinear emissions, including those not singular as Zi — > 1. Because it contains the terms 
which dominate at large p±, as well as those that dominate at small one may hope 
that this approximation will be good in the intermediate range of p± as well. 
In analogy to Ref. [16| we define the ratio 



^sing ^ {da/dp\/dyH) ("singular" only) ^^^^ 
[da / dp\/dyH) (exact) 



We plot the ratio R^™^ in Fig. |10| for two choices of the Higgs boson mass. We also plot 
the similarly-defined ratio i?^"^^, where we have calculated the cross section with soft-plus- 
virtual terms only, by keeping only those terms in eqs. (|3l7|) , (|3^ ), d^), and (|3^ that 



are multiplied by 5{Q'^) or a +function singularity, and setting = in all factors that 
are multiplied by a +function. Note that since we write the +functions in terms of Zi rather 
than Q^, our definition differs from [|^] by terms that are nonsingular as — > 1. From the 
figure, we see that the pure S+V approximation does an excellent job at p± = 200 GeV for 
TJiH = 120 GeV, overestimating the exact NLO cross section by 2.5%, while the "singular" 
approximation also does well, underestimating by 6.7%. The surprisingly good agreement 
of the S-fV approximation is somewhat fortuitous; keeping some non-leading dependence 
on Zi typically changes the cross section by about 5 — 10% at p± = 200 GeV. At medium to 
lower p±, the "singular" approximation does better than S+V. For instance, at p± = 100 
GeV the "singular" approximation underestimates by 4.7%, while the S+V approximation 
overestimates by 11.6%. This suggests an explanation for the observation in Ref. |l^ that 
the S+V approximation does worse for larger uih- Although increasing mn does restrict 
the phase space for real emission, it also makes the logarithm log niH /p± larger for a given 
p±. Thus, the range in p± for which the small-pj_ logarithms dominate increases with 
rriH- To illustrate this effect, we have also plotted the same ratios for mn = 200 GeV. 
For this larger Higgs mass, the "singular" approximation improves, whereas the pure S+V 
approximation gets worse. 



6. Conclusions 

In this paper we have presented a NLO calculation of the Higgs boson transverse momenta 
and rapidity distributions, for nonzero Higgs p±. We have included all of the analytic 
formulae, so that it is possible to directly investigate them in various kinematic limits. In 
particular, we have checked the small-p^ limit, and directly verified the formulae predicted 
from the resummation of large logarithms of Inmn /p±, including the B^'^^ coefficient. 
We have also isolated the most important pieces of the cross section into the "singular" 
contributions, which include the soft+virtual contributions, the remaining collinear contri- 
butions, and most of the small-p_|_ contributions. This may be useful as a starting point 
for a combined soft+virtual and small-px resummation. 

The numerical size of the NLO corrections are large, with a iT-factor rising slowly from 
1.6 to 1.8 as the p± increases from 30 GeV to 200 GeV. The scale dependence is reduced at 
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NLO, compared to LO, with a variation on the order of 20% for 0.5m± < /i^ = fip < 2m±. 
We expect that the theoretical uncertainty from uncalculated higher orders (as well as due 
to the nit ^ oo limit) is larger than that due to uncertainties in the PDFs, although a 
detailed error analysis on the effect of these uncertainties on this particular process (as 
outlined in |3^]) would be necessary to make this definitive. 
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A. Nonsingular Real Contributions 

Here we give the "nonsingular" contributions coming from the real H + 2 parton cross 
section. These terms are finite as e ^ and as 0, so that they can be written in 

4 dimensions. We have also removed any naive singularities in p^, before the phase space 
integration of eq. ( |2.3| ). Of course, this separation into "nonsingular" and "singular" terms 
is not unambiguous, but the sum is definitely unique. We find the separation useful, since 
the formulae for the "nonsingular" terms are very long, whereas the dominant contribution 
to the cross section comes from the "singular" terms given in the main text. 

A.l gg ^ H + X terms 

We give the result as a sum over terms coming from specific color-ordered helicity ampli- 
tudes (after subtracting from each the "singular" terms that contribute to Gj-^^'*^^), plus a 
correction term arising from the 0{e) correction to the matrix element times a collinear sin- 
gularity, which is necessary to obtain the correct CDR result. That is, for the gluon-gluon 
"nonsingular" contribution we write 



Gf^'''"'^ - -2^5 g2|^c (^0+^(1234) + ^(2341) + ^(3412) 



+^(4123) + ^(1324) + ^(2413) + ^(3241) + ^ej (A.l) 
+n/Cp(5i(+_) +n/iVe(52(+_) +52(++)) \ ■ 



For definiteness, Aq comes from Hgggg squared amplitudes with helicities (-1- -|- -|— 1-), ( 

), or {+ -\ ), where all helicities are taken to be for out-going momenta. The ^(ijfci) 

come from the Hgggg color-ordered amplitude m{i—,j-\-,k+,l+) squared, as well as all 
other squared amplitudes that are identical by color-reversal, helicity-reversal, or exchange 
of the final-state gluons. For example, ^(3412) comes from m(3— , 4-|-, 1-|-, 2-|-) plus all 
identical color-ordered amplitudes, where 1 and 2 refer to the incoming gluons. (We use the 
cyclical symmetry of the color-ordered amplitudes to always put the negative helicity gluon 
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in the first position.) The terms, ) and come from Hggqq squared amphtudes 

with unhke-hehcity gluons and hke-hehcity gluons, respectively. Finally, comes from 
the C(e) corrections to the squared Hgggg matrix element that are not proportional to 
(5((5^) and cannot be identified as a contribution to the C|j, and therefore have not been 
included in G^gg \ eq. ( 3.17| ). The Hggqq squared matrix element also has 0{e) corrections, 
but they have all been included in Gf^\ eq. ( ^Tzl) . 
For later use in these formulae, we define 
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where 

^ = s + - Q2 , 
B = -4m^s)^ 
We obtain from the Hgggg squared matrix elements: 
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+mjjQ^{134i + llO^i) + 4Q^ + 52mjiQ^ + 20m%Q'^ - 22m%Q^ 

-4(4123) = -4(3412) with |(t,a) ^ {u,b)^ , 
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We obtain from the Hggqq squared matrix elements: 
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A. 2 gq ^ H + X terms 

We write the gluon-quark "nonsingular" contribution as 

(-»{2R,ns) 1 
99 



s2p^Q2 



{— ) 



(A.17) 



+NcCf\C2(+-) + C'2{-+) + C2{++) + C'2{ ) + C2e 



where the terms, Cj(±-|_) , come from the Hggqq squared amphtudes with hehcity configura- 
tions {q— ■,(!+■, 91 ±, 32=1=); plus the hehcity-flipped configuration that contributes identically. 
To be precise, the helicities are written for all final-state momenta, before crossing the par- 
ticles gi and q into the initial-state to obtain the scattering giq — > g2qH. The term 
comes from the 0{e) corrections to the squared Hggqq matrix elements that have not been 
included in cff, eq. (Isl^l) . Of course, Gg^^''^''^ is obtained from these expressions by 
(t, a) <-> (n, 6). 
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We obtain: 



C 



a 



2(-+) 



c 



C2(+_) = , 



l2 X\3 



Pi 



2Qi 



(|(Q2_£)3 + Q2(g2_^)3) 



KQ^-t? + g^(Q^-u)^) (-2 + 3^) 



Ss^'pI spIQ^A^ 



2u 

Ml, 



uB 



+^spi{Q'-tr 



+^sp'i{Q^-u] 
§pl{Q'-i) f^ 



(g2 - n) u 



± 

2 /n4 



(Q2 -t) it 
2su{s + £) - Q'^{4m%Q'^ - QH - mjfU 



-2ui) 



-1 + 2^^ 



ts^ - 2uts + 2Q^u{Q^ - t) 



spjiQ'-u) 
u{Q^ - i) 

^sp']_Q'^mjj{i + u) 2sp\Q'^m'\j ^ s^p\Cfrn\j 



-1 + 2 



2u 



-2((5^ + m]j)su + 2n?HS^ + uihS 
+m%Q^(2{s - Q2) + smjj ~ uj + 5Q'^s{s - 



(A.18) 
(A.19) 



(A.20) 



(A.21) 



-L2b 



{A - mjjYsr - 2Q^mijutA - Q^rn\j{Q - 0« 



{s + Q^- ml) - mlf + Q^A^'^ ~ 4sQ^A{A - m]j) 



(A.22) 



C2(++) = T^SpiA\l - Za){L2a - L^a) + 



spl{ml-i)A^{l-Zi) 



2u 



{Lib — L2b) 



2u 



L 



lb 



{mjj - u 

2u{s - Q^f + Am\{m\ - t)A - 2m%{Q'^ - i) - m% 



IspiiQ'-ty 



(g2 -u) u 



Q2 Q4n 
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Ci(— ) 



+ 



Q2 



2u(Q2 _ 



(Q2 -t) U 



2u(Q2 _ i) 



3spl{s + Q') - uQ\Q' -t) + SQ'mjjiQ' - u 



+st{u + s) 



-1 + 2 



2 ^2 r)2 



sp±mjjQ 
2^ 



2(s - Q'y - 2mjj{s - mjj) - u{Q'^ - u) - AmjjQ^ 



s^p\P 
u 



L 



la 



2u 

sp\Q'^{m]j-if 



-Lib , (A.23) 



u 



C2(__) = 



^jiy^' - mjjf{u + t) + 2Q^m% + Q^m%{^i - u) 
{-Ku + i) + 2mlQ^ + Q\i - u)) , 

+ 3Li5) + \spif{l - Za){L2a - ^la) 



(A.24) 



2^2 



(m|f - u) 



2(sj)im2,g2)2 



sp\Pm?jjQ'^ 8{spj_mffQ 



S^2 ^2 n2\2 



u4 

2 ^2 r)2\2 



U4 



C2e = 



4:{spj_mjfQ^) 



(A.25) 
(A.26) 



A. 3 QiQi H + X terms 

We write the quark-antiquark (same flavor) "nonsingular" contribution as 



G. 



(2R,ns) 



+nfClEi + ClE2 + ^EA , 



(A.27) 



where the terms, D^^^ ■) and D^^^_^^^ , come from the Hggqq squared amphtudes with unhkc- 

helicity gluons and like-helicity gluons, respectively. The Ei terms arise from the s-channel 
diagrams, the E2 terms come from the t- and ^-channel diagrams, and the E^ terms 
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come from the cross terms. In the case of quark- antiquark and quark-quark scattering 
contributions, all 0{e) corrections have been included in the formulae in the main text for 
off, eqs. (PI) and (p|). 

We obtain from the Hggqq squared matrix elements: 



n 



2(+-) 



D. 



2(++) 



t 



+ 



t u 



(A.28) 



a^2 ;2 



+ 



t U 



(A.29) 



A 



2b — ^Ib 



Lib) + 



A 



{L2a — Lia) 



sp\rnyQ^{sp\ + ut) 2s^p\rn\jQ'^ sp]_mjjQ'^ {2sp]_ - ui) 
H ^ — L3 ' 



t2 



AB 



+ 



t ^ u 



(A.30) 



i-t^lh - ^2b) H ;r-; [^la — ^2a) 



2A 



2A 



--Sp_L« [Lia + Lib) H ^4 ^2 



f Sp\v?[u + t) , 2 2 ^2\ 



2 _2 ^2\2 



+ 



B 



+sp\m\jQ^ 



t ^ u 



1 



3^ 



AB^ B^ 



2B 



4^3 



(A.31) 



We obtain from the Hqqqq squared matrix elements: 



771 ° '2 2 ,o2 ^ - 2 2 ,o2 
^1 = P±Q - -^sp^m^Q , 



E2 
E3 



spiQ\Q'-pi) 



{it/zaf + {u/zhf)+2sV^Q' 



-2spl ({u + i- 2Q'^ f - 2spl) Inpi/Ql 



spiQ\2Qi + Q') 



(A.32) 
(A.33) 

(A.34) 



A. 4 qiqj —^H + X terms 

The quark-antiquark (different flavors) "nonsingular" contribution comes from Hqiqiqjqj 
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squared amplitudes and can be written 



G. 



(2R,ns) 



CpE2 



where E2 was given in eq. ( A.33| ). The quark-quark (different flavors) "nonsingular" con- 
tribution is also given by this same expression, eq. ( |A.35D . 

A. 5 qiQi —^H + X terms 

The quark-quark (same flavor) "nonsingular" contribution comes from Hqiqiqiqi squared 
amplitudes and can be written 



^(2R,ns) 



1 



where E2 was given in eq. ( A.33| ), and the subleading term in \/Nc is 



Q] 



(A.35) 



B. Phase space for integrals over momentum fractions 

After integration over the angular variables of the two emitted partons, we are left with 
the integrals over Q < Xa,Xh < ^ with the restriction to > 0. A parametrization of this 



phase space was given in Ref. |40]. Using a change of variables, we can write the phase 
space as a sum of two double integrals, where the inner integral is over the momentum 
fraction Za or Zfe, defined in eq. ( |3.8| ). Specifically, for the case where the inner integral is 
over Za, we write 



1 

dXa I dXbQ{Q'^)F{Xa,Xb) 

Jo 



1-5 



dzi 



dZa 



ZaZ'b 



+ 



1-5 



dz' 



dZa 



Z„ ./x.0/4 Za 



p\ 



(1 - z'aYza 



(B.l) 



F{Xa,Xb) 



with 



In the first double integral 



m_|_e 



VH 



(B.2) 



^(l + Aft) 
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2 



2 

m 



^(1-4)(1 + Afc) 



u = mjj + 



2 

mi 



= ^(l-^„)(l-4)(l + A6) , (B.3) 

where 



ml (1 - 



(B.4) 



In the second double integral 



^0 



Xb = Xb{l + Xa) 



pi 



u = mjj 



m 



2 



For very large | yn \ the limits on the integrals may be further restricted by the requirements 

Xa-} Xb ^ !• 

The two double integrals become particularly simple in the limit p±/m± 0, and they 
have simple physical interpretations. The first double integral becomes two independent 
convolutions, corresponding to independent emissions off the a and the b partons. The 
second double integral becomes a double convolution, corresponding to two successive 
emissions off the a parton. We can alternatively write the phase space with the inner 
integration over Zb by replacing a <^ b, t ^ u everywhere. The choice of the inner integral 
over Za or Zb is a matter of convenience. We have found it most natural (and efficient in the 
small-pj_ limit) to compute the "singular" terms that are written specifically as functions 



or distributions in Za{b) in section |3.4| using the phase space with the inner integral over 



Za(b) ■ The other terms are computed with the phase space symmetrized over a <->■ b. 
C. Details of the small-p^ limit 



As discussed in section |5|, the resummed formula, eq. (5J), can be expanded as a power 



series in as, yielding 
da 



where the coefficients are 



s p\ 



■ 2 2m-l 

EE 

.m=l n=0 



V27r 



m 



H 



iCn In 2 

PI 



+ 0{al) 



(C.l) 



iCi 
2C2 



A(i)/,(xO)/,(a;°), 



1 

2 L 
3 
2 



^(1) 



2C1 



^(1)1 2 _ ^(2) _ ^^^(1) _ 3^^^(1) 1^ i4 



H 



fgK)fgi^t) 



+{a ^ b} 



2C0 



Pg^ ° fi) i<) (Pg, ° f,) (4) + (Pg^ « ° Z.' ) i4)f9i4 



In 



+ 



3^oln^-i?«ln^ 



+ 



9J 

C(3) [^«]% V) + ^/3oB«ln5 



.(1) 



9^ 



m 



H 



fgi4)fg{4) 



+{a ^ 6} . 



(C.2) 



In these equations the indices i and j are implicitly summed over all massless partons 
{g,qf,qf), and the parton density functions are always evaluated at the factorization scale 
Hf', that is fi{x) = fi{x, ^p)- The parameters, A^^'"^"*, 2?^^'^), and the functions clj\z) were 
given m eq. Q. The function p(^\z) is the NLO splitting function. The expressions for 
the mCn are precisely the same as those given in Ref. ||3^, except for the additional factor 
of 3 in front of each occurrence of ln(;U^/m|^). This factor, of course, is due to the fact 
that the LO p± spectrum for gg ^ H + X with > begins at ©(a^). 

We have explicitly checked that our NLO calculation agrees with this asymptotic ex- 
pression in the limit p±_ <C mn- We give some of the details of how this check was performed 
here. 
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The coefficients in eq. ( p.2[) can be separated into contributions: 



(C.3) 



where the superscripts label which type of partons come from the a and b hadrons, and 
the label q actually implies a sum over all light quarks and antiquarks. The quark-quark 
contributions, mCn''\ are particularly simple, so we begin the discussion with them. 

First we note that Qqq is not singular as p± 0, so that the LO cross section does not 
contribute to the small p± limit; i.e., iCn'^^ = 0. Similarly, all terms at NLO which are 
proportional to 6{Q'^) do not contribute. In fact, it is easy to see that the "singular" terms 
that contribute are independent of quark flavors and of whether it is quark-quark or quark- 
antiquark scattering. Finally, we note that none of the "nonsingular" quark-quark terms 
contribute in the small-p^ limit. The general argument is as follows. Since the Q^^^'^^^ 
terms are defined to have no explicit singularities as p± — > 0, they can only contribute at 
if they obtain the singularity through integration over the phase space near = 0. 
Effectively, this implies that "nonsingular" terms can only contribute if they are 0{l/p\) 
when is treated as ^^'^ either n or t is treated 0{p\). 

Thus we get 



1 



-P,,(z,)lng^ + q,(z,) 
(s-Q')' + (^i + t-2Q2) 



^In^ 

pI Ql 



+ 



{t, a) ^ (n, b) 



(C.4) 



Inserting this into eqs. ( |2.4D and (p.Sp, this is easily integrated using the phase space 
parametrization of eq. (B.l). Since there are no singular terms in the integrands as ^ ^ 1 
or Za^b ~^ 1; the p± <^ niH limit is trivially taken, with only the first integral of eq. (B.l) 
contributing. The result is 



2C. 



(qq) 



2C' 



(qq) 
2 







(qq) 



2 Co 



(qq) 



Pm ° fq) {<) [Pgq ° fq) (4) + {a ^ b} , 

[Pm o fq) (x°) (P,, o /,) (4) In £f + o /,) (x°) (p,, o /,) (x°) 

+{a ^ b} . (C.5) 



The gluon-quark contributions are obtained similarly, with two added complications. 
The first is that the "nonsingular" terms now contribute to the cross section. The simplest 
way to handle these terms is to note that they are singular as u ~ Oi^p"^) — > but not 
as t ~ ©(p^) ^ 0. Thus, we can use the phase space of eq. ( |B.l| ) with the replacement 
a ^ b, and only the second integral contributes in the small-p^ limit. Since there are no 
subtleties in this integral, one can directly take the limit p±/mH everywhere, leaving 
a double convolution, which can easily be converted to a single convolution. 
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The second complication is that the integrands have singular behavior as both z'^i^ ^ 1 
and as Za^b 1- The first of these can be handled by the following identity as p±/mH — > 0: 



In'^fl - 
i — z'- 





"ln^-(l-z9' 


/ dz'i F{zi) 

JxO 





1 2 \ fc+l 

1 , mTr 
--In — - 

2 -2 



PI 



(C.6) 



where we recall that b = p^/{mx_ +p±). For the second case of inner integrals of the form 
we require the following identities as p±/mH — > 0: 
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vr 



21 



6 



(C.7) 



where A and A are 0{p\/m?^. 

The gluon-gluon contributions have only one additional complication beyond these. 
On integrating lines 5-6 and lines 9-10 of eq. ( |3.17| ), one encounters double convolutions, 
which can be transformed into single convolutions. However, extra care must be taken 
with the endpoints of the integrals as p±/mH 0. In the process we require the following 
identities in this limit: 

1 , , /1-z'- 



l-<5 

,0 



dzlFizDj 



In' 



1 



2C(3)F(1) , 
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dz' F(z') 
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9 2 

ln2^ +2C3 

12 



F(l). (C. 



We have verified that the coefficients mCn'^^ and mCn^\ derived from our NLO cross 



section agree with that of eq. ( p.2| ). As an example, we display separately the "singular" 
and "nonsingular" contributions to 2Cq^ . The singular contribution gives 



2C0 



P99 ° fa) {<) (P99 ° fa) + [Pga ° Paa ° fa) iOfaiO 



In 
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H 
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where 
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The double convolutions have been transformed into single convolutions in the second form 
foi Diiiz). 

We can write the "nonsingular" contribution as 
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The convolution function is a sum over terms coming from each of the squared hehcity 
amphtudes: 



+^[A(1324)+A(2413)](^) + ^A(324l)(^) + ^A, 



(C.ll) 
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where the labels correspond to the terms in eq. ( [A.l| ). We obtain 
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In 2 , 

(C.13) 
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Note that the 0(e) term contributes and is crucial to obtaining the correct small-p^ 
hmit. 

Combining sC^^'^^^ + 20^'"'^ and comparing with the coefficient in eq ( |C.2| ) obtained 
from the resummation, we recognize the combination 



b(^)6{1 -z) + P^f{z) = D^;Jiz) + D^^;Hz) 



(l-z)_ 



1 67 5 vr 
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This agrees with the result for i?^^) from Ref. |2C] and the result for (z), given in 
Ref. ||. 
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Figure 1: The Higgs p± spectrum for tjh — calculated at LO, with the exact to* dependence 
(dashes), and in the rrit ~ oo effective theory (solid). 
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Figure 2: The Higgs p± spectrum at NLO in the rat — oo effective theory for rapidities uh = 
(solid), yn = I (dashes), and yn = 2 (dotdash). 
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Figure 3: The Higgs p± spectrum for = at NLO in the nit — oo effective theory, separated 
into its initial-state partonic subcomponents: gluon-giuon (dashes), (anti)quark-giuon (dotdash), 
and total (solid). The (anti)quark-(anti)quark contribution is very small and negative over the most 
of the range of this plot. 
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Figure 4: The Higgs p± spectrum for yn = a.t LO (dashes) and NLO (solid), with the renor- 
malization and factorization scales set equal fiR — fip and varied in the range 0.5m± to 2m±. 
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Figure 5: The if-factor for the Higgs pj. spectrum at uh = 0. The scales are set to ^i?^ = = wj.- 
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Figure 6: The Higgs cross section for yn = Q and p± = 50 GeV as a function of varying renor- 
malization fin and factorization fip scales, normalized to the cross section at /in = hf = m±. The 
scales are varied together at LO (dotdash) and NLO (solid) and also inversely at LO (dots) and 
NLO (dashes). 
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Figure 7: The rapidity dependence of the Higgs cross section for two different transverse momenta. 
The curves are for pj_ = 50 GeV at NLO (sohd) and LO (dotdash) and for p± = 100 GeV at NLO 
(dashes) and LO (dots). 
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Figure 8: The i^T- factor as a function of rapidity for p± = 50 GeV. The scales are set to /i/j 
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Figure 9: The Higgs p± spectrum (solid curves) compared to the small-p^ limit formulae Eq. (5.8) 
(dotted curves) at both LO and NLO. All curves are at calculated at y^f = in the mt — oo effective 
theory. 
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Figure 10: The ratio of the NLO p± spectra, calculated with cither "singular" terms only or 
soft+virtual terms only, divided by the full NLO cross section. The curves are for mn — 120 GeV 
(solid) and niH = 200 GeV (dashes). All cross sections are for y// = in the mt = oo effective 
theory. 



-40- 



